Perfect quantum transport in arbitrary spin networks 
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Spin chains have been proposed as wires to transport information between distributed registers 
in a quantum information processor. Unfortunately, the challenges in manufacturing linear chains 
with engineered couplings has hindered experimental implementations. Here we present strategies 
to achieve perfect quantum information transport in arbitrary spin networks. Our proposal is based 
on the weak coupling limit for pure state transport, where information is transferred between two 
end-spins that are only weakly coupled to the rest of the network. This regime allows disregarding 
the complex, internal dynamics of the bulk network and relying on virtual transitions or on the 
coupling to a single bulk eigenmode. We further introduce control methods capable of tuning the 
transport process and achieve perfect fidelity with limited resources, involving only manipulation of 
the end-qubits. These strategies could be thus applied not only to engineered systems with relaxed 
fabrication precision, but also to naturally occurring networks; specifically, we discuss the practical 
implementation of quantum state transfer between two separated nitrogen vacancy (NV) centers 
through a network of nitrogen substitutional impurities. 

PACS numbers: 03.67.Ac, 03.67.Hk 



Transport of quantum information between distant 
qubits is an essential task for quantum communication [1] 
and quantum computation [2]. Linear spin chains have 
been proposed [3] as quantum wires to connect distant 
computational units of a distributed quantum processor. 
This architecture would overcome the lack of local ad- 
dressability of naturally occurring spin networks by sep- 
arating in space the computational qubit registers while 
relying on free evolution of the spin wires to transmit 
information among them. Engineering the coupling be- 
tween spins can improve the transport fidelity [4], even 
allowing for perfect quantum state transport (QST), but 
it is difhcult to achieve in experimental systems. Remark- 
able work [5, 6] found relaxed coupling engineering re- 
quirements - however, even these proposals still required 
linear chains with nearest-neighbor couplings [7, 8] or 
networks will all equal couplings [9] . These requirements 
remain too restrictive to allow an experimental imple- 
mentation, since manufacturing highly regular networks 
is challenging with current technology [10-13]. 

Here we describe strategies for achieving QST between 
separated "end" -spins in an arbitrary network topology. 
We employ the weak-coupling regime [5-9, 14], where 
the end-spins are engineered to be weakly coupled to the 
bulk of the network. We describe the transport dynam- 
ics via a perturbative approach, identifying two different 
regimes. Perfect transport can be achieved by setting 
the end-spins far off-resonance from the rest of the net- 
work - transport is then driven by a second-order pro- 
cess and hence is slow, but it requires no active control. 
Faster transport is reached by bringing the end spins in 
resonance with a mode of the bulk of the network, effec- 
tively creating a A- type network [15], whose dynamics we 
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FIG. 1. Example of spin network, consisting of NV centers 
(blue spheres) and PI centers (red) in a diamond lattice. The 
network is given by randomly positioned PI centers in the 
diamond lattice with a concentration of 0.2 ppm and a 5% 
conversion efficiency to NV. The proposed strategies enable 
perfect quantum state transfer between the two NV spins in 
this naturally occurring topology of PI centers [20, 21]. 



characterize completely. We further introduce a simple 
control sequence that ensures perfect QST by properly 
balancing the coupling of the end-spins to the common 
bulk mode, thus allowing perfect and fast state trans- 
fer. Finally, we investigate the scaling of QST in various 
types of networks and discuss practical implementations 
for QST between separated nitrogen vacancy (NV) cen- 
ters in diamond [16, 17] via randomly positioned elec- 
tronic Nitrogen impurities [18-21]. 

Spin Network - The system (Fig. 1) is an A^-spin net- 
work, whose nodes represent spins- ^ and whose edges Hij 
are the Hamiltonian coupling spins i and j. We consider 
the isotropic XY Hamiltonian, Hij = {Sf S~ + S~ Sj'), 
with sf = ^{Sj ± iS^), which has been widely stud- 
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FIG. 2. (a) On-resonance balancing sequence applied to the 
end spins for perfect quantum transport. Yellow blocks indi- 
cate microwave irradiation that brings the end spins on res- 
onance with the bulk network. The resulting A-network is 
in general unbalanced, but appropriately placed 7r-pulses (or- 
ange) on the end-spin with higher mode overlap (here spin 1) 
can balance the overlaps 0{i,jv) [25]. (b) The vr-pulses invert 
the sign of spin-1 coupling to the bulk mode in the toggling 
frame such that Oi,jv become equal on average. For the net- 
work of Fig. 1, |Oi|>|Ojv| andr=l/2(l-H|Oiv/Oi|) = 0.5501. 
The sequence is symmetrized [26, 27] and repeated for L cy- 
cles. 



sibly degenerate) mode. The system reduces to a A- 
network [15], where coupled- mode theory ensures perfect 
transport if both ends have equal overlaps with the bulk 
mode [29, 30], a condition that we will show can be en- 
gineered by the control sequence in Fig. 2. 

To make more rigorous our intuition of the weak 
regime, we describe transport via a perturbative treat- 
ment. For convenience we consider normalized matrices, 
\\Hb\\ = ||i?e|| =1, and introduce the network adjacency 
matrix, A — f3A-^ + eA^, which describes the coupling 
networks of the system Hamiltonian H ~ I3Hb + ei/g- 
Transport in the single excitation subspace is fully de- 
scribed by A [4], thus the fidelity can be written as 
F{t) = I (7V| e"*"^* |1) 1^, where the vectors \j) now repre- 
sent the node basis in the N x N network space. We use 
a Schrieffer- Wolff transformation [31, 32] and its trunca- 
tion to first order in e/(3 to define an effective adjacency 
matrix. A' — e^Ae~^ « A+ ^[S,A], which drives the 
evolution. Setting S so that [^^,5*] — ^A'^, we have 

A' sa (3A^ + ^A^, where A^ = [S, A"] can be evaluated 
explicitly. 



ied for quantum transport [3, 15, 22, 23]. We further 
assume that two nodes, labeled 1 and iV, can be par- 
tially controlled and read out, independently from the 
hulk of the network: we will consider QST between these 
end spins. For perfect transport, an excitation created 
at the location of spin 1 should be transmitted without 
distortion to the position of spin N upon evolution un- 
der H . We characterize the efficiency of transport by the 
fidehty, F{t) ~ \ {N\e'~^^*- |1) |^, where \j) represents a 
single excitation |1) at spin j, while all other spins are in 
the ground state |0). 

Weak- coupling regime - While optimal fidelity has 
been obtained for particular, engineered networks 
(mainly ID, nearest-neighbor chains), here we consider a 
completely arbitrary bulk network, Hb- To ensure per- 
fect transport, we work in the weak-coupling regime for 
the end-spin coupling H^. By engineering appropriate 
weights e,/3, we thus impose e||^^e|| ^ /3||-H^b||, where 
II • II is a suitable matrix norm [24]. 

Intuitively, we expect the weak coupling regime to 
achieve perfect transfer since it imposes two rates to the 
spin dynamics: the bulk spins evolve on a "fast" time 
scale while the end-spin dynamics is "slow" . The end 
spins inject information into the bulk, which evolves so 
quickly that information spreads everywhere at a rate 
much faster than new information is fed in, allowing an 
adiabatic elimination of the information quantum walk 
in the bulk network [15, 28]. Although high fidelity can 
be reached, this off-resonance transport is very slow. 

A different strategy, and a faster rate, for informa- 
tion transport is achieved by bringing the end-spins on 
resonance with an eigenmode of the bulk - the weak 
coupling ensuring greater overlap with a single (pos- 



Off-resonance QST - Consider first the case where 
the eigenvalues {E^} of A^ are non-degenerate, except 
for — Efj = (associated with the end-spin subspace) . 
We can fix the energy eigenbasis {[wfc)} of A^ by setting 
|ui) = |1) and \vn) = \N). In this basis the structure 
of the matrix A^ is preserved, non-zero terms connecting 
only the ends to the bulk, A^ - = Q ior i ^ l,N. A 
general element of A^ can be written as. 



,2 / -1 

AS _ 1_ AC AC \ 



Ef-E^ Ef-E^ 



Given the form of A^ ^ we have, for {«, j} 7^ 1, iV, 



A% = ^(A.iAy + A^atA^,) f ^ + j ■ (1) 



Also, setting fc 7^ 1, iV and {C, ^} e {1, N} we have 



2£l 
/3 



(2) 



while if A| ^ — elements between the end and bulk 
are zero, — 0. Hence A^ can be partitioned into 
a term with support only in the bulk subspace (Eq. 1) 
and a second term with support only in the end-spin 
subspace (Eq. 2), while there is no end-bulk coupling. 
To first order approximation, it is only the latter term 
that drives QST - the transport happens via a direct 
coupling between the end-nodes. Since this effective cou- 
pling is mediated by the bulk via virtual transitions, its 
rate is proportional to e^//3. The fidelity of transport is 
determined entirely by the effective detuning of the two 
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FIG. 3. (a) Off-resonance transport fidelity for tlie networlc in 
Fig. 1 for 7 = 25. Perfect transport occurs, but on a slow time 
scale, an order of magnitude longer than the on-resonance 
balanced case (Fig. 4). (b) Increasing 7 improves transport 
fidelity (red circles) but also increases the time required for 
perfect transport (blue diamonds). 
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FIG. 4. (a) On-resonance transport fidelity for the network 
of Fig. 1, with (red) and without (blue, dotted) balancing, 
with 7=1 and L = 20. Almost perfect transport occurs 
in the balanced case, obtained by the control sequence in 
Fig. 2. (b) Increasing the number of cycles L improves the 
Trotter approximation yielding enhanced fidelity. The sym- 
metrized sequence performs better than the sequence without 
symmetrization. 



end-spins, a = {Af-^^ — A|,^)/2: 

m = i^AS^yla^ ^i^' {t^{Af^Y+a^^ (3) 

If we can modify the end-spin Hamiltonian by adding a 
term iJoff — —2{ujxSl + ujnS^), such that Af-^^ -|- = 
■^fiN + ^JV, perfect quantum transport is ensured at 
trn — ■''■/(2^fjv)- This energy shift could be obtained by 
locally tuning the magnetic field or by applying local AC 
driving, ensuring the desired energy in the rotating frame 
(similar to the Hartman-Hahn scheme [33]). Transport 
fidelity also depends on the goodness of the first order ap- 
proximation, increasing with /3/e as shown in Fig. 3(b) 
at the cost of longer transport times. 

On-resonance QST - Transport can be made faster if 
the end spins are on resonance with one non-degenerate 
mode of the bulk \vd) (we will consider the degenerate 
case in [34]). Resonance happens if the corresponding 
eigenvalue = or it can be enforced by adding an en- 
ergy shift to the end spins to set A^^ ~ A%j^ — f3/eE^ . 
Transport then occurs at a rate proportional to e, as 
driven by the adjacency matrix A'^ ^ the projection of A^ 
in the degenerate subspace: 

A'' = (1] A'^ \v,) \l){v,\ + {N\ A^ \va) \N){v,\ + h.c. 

We note that the goodness of this approximation depends 
on the gap between the selected resonance mode and the 
other bulk modes. In the node basis, A'^ forms an ef- 
fective A-network [15], coupling the end-spins with each 
spin of the bulk: 

A'' = J2{S,,\l){j\+S^,\N){j\+h.c.), (4) 

j 

where (5(i,jv)j — {''^d\j) {{^, N)\ A^ \vd) ■ Note that im- 
portantly we have Sij/Sj^j = est., Vj. Transport in 
such A- networks occurs at only four frequencies [34], 



F{t) = Wo + J2m=i cos(/,„t) with, 

/i,2 - 2^S^T^S^^ ; /3,4 = V2(^2tA2), (5) 
where, 

3 

= I](^ij-4iv - 5jN5ikf (6) 
j<k 

3 

Physically, S ^ ||^^|| sets the energy scale of the resonant 
mode, while A quantifies the relative detuning between 
different A-paths [15, 34]. The coefficients Wi are found 

to be Wo = -W3 = -Wi = 4(g4,A-*) ' wi = W2 = wo/2, 
giving the analytical expression for QST in a A-network, 

Fit) ^ sm^ (y ^7 IV 

(7) 

Perfect QST requires A = and 6 = S. The first condi- 
tion entails Sij/SjN = est. Vj, which is always satisfied by 
A'^ if the resonant bulk mode is non-degenerate. On the 
other hand, S = S requires a balanced network, Sij — Sjn 
Vj. For the reduced adjacency matrix A'^ this condi- 
tion is satisfied when both end-spins have equal overlap 
with the resonant eigenmode, (l|A'^|z;d) = (-^1 I'^d): 
and we show below that this can be always arranged 
for a non-degenerate mode by a simple control sequence 
(see Fig. 2). In the balanced case the fidelity simplifies 
to F{t) = sm{St/V2)^, which leads to perfect QST at 
tm = 7r/(\/2S'), as if it were a 3-spin chain [4]. 

Perfect QST by on-resonance balancing - Unfortu- 
nately the overlaps of the two end spins with the on- 
resonance mode, 0(i_7v) = ii^: ^)\^'^\vd) , are in gen- 
eral unequal, and the A-network unbalanced. Still, it 
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is possible to balance the network by a simple control 
sequence (Fig. 2). Assume for example Oi > On', we 
partition into effective adjacency matrices with cou- 
plings only to spins 1 and N, A'^ — Af + A"^. A ro- 
tation e""^^'^ on spin 1 produces A''- — —Af + Af^. 



iAtt 



with 



Then the evolution, e 
Af^ = A'^r + — r), is balanced on average during 
the period t ifr = ^(1 -|- \0n/0\\). The approximation 
improves if one uses L cycles of the control sequence with 
shorter time intervals (as in a Trotter expansion [27, 35]) 
and appropriately symmetrizes it (see Fig. 2) to achieve 
an error 0{r'^{l — r)t'^/ L^). Fig. 4(a) shows the effect 
of enhanced, almost perfect, fidelity upon balancing the 
network of Fig. 1. 

Transport time and control requirements - We now 
consider the scaling of the weak coupling parameter 
J — J required for the validity of the perturbative ap- 
proximation. We can fix /3 = \\Hb\\ and consider nor- 
malized matrices. Then /3 scales as i/^^V^^^ypV^^^JTS 
for a random network where all the couplings are uni- 
formly distributed. For the more realistic case where the 
coupling strength decreases with distance, the scaling is 
less favorable, e.g. for a random dipolar coupled network 
f3 scales as y/2/3{N — 3)/dP, where d is the average lat- 
tice constant [34]. Similar scaling occurs for regular spin 
networks, for example those corresponding to crystal lat- 
tices [34]. In general the ratio e//3 decreases with the size 
of the network (usually as 0{\/N)), averting the need to 
reduce the end-couplings by engineering e. This is evi- 
dent in Fig. 4, where e = 1 is sufficient to drive perfect 
quantum transfer. 

In the case of on-resonance balancing the time at 
which perfect QST is achieved is tm = JT^/iV^S), where 
S = min{Oi,OAr}, scales linearly with 7. The time 
is shorter the more symmetrical the end-spins are with 
respect to the resonant mode [34], since then |Oi| ~ 
I Oat I . For the off-resonance case the time is t^ — 



{7rf3/2e^)[Ef/{Al,Al,)], 



where Elr = 



fc |} is the 



minimum eigenvalue of the bulk. In general this sec- 
ond order transport process takes an order of magnitude 
longer time than the on-resonance case (Fig. 3). 

Finally, let us estimate the resources needed to impose 
the end-spin energy shifts as required for perfect QST, for 
example by a continuous irradiation during (Fig. 2). 
In the on-resonance case, the end-spin energy should be 
set equal to a bulk mode, . Selecting the highest 
bulk eigenmode, which has usually the largest gap to the 
other modes, E^ scales as 0{N) for random networks, 
but it is constant, ~ 1.6/d^, for dipolar coupled networks 
[34]. Off-resonance transport requires instead a shift of 
- ^NN^ where Af^j^j^ = {Alj^\if /j^Ef , where 
Ef is the minimum bulk eigenmode. Hence the control 
required in this case is about an order of magnitude lower 
than in the resonant case. 



An experimentally important quantum computing ar- 
chitecture is that of nitrogen vacancy (NV) centers in di- 
amond. Fig. 5 shows the scaling of the transport time be- 




o 



"16 20 24 28 

Distance between end-spins (nm) 

FIG. 5. Time for optimal transport in dipolar coupled spin 
networks of PI centers in diamond. The optimal time was 
calculated from the average over 5000 random lattice realiza- 
tions of PI centers of density lOppm; we considered transport 
between two NV centers located at increasing distance. The 
red line (circles) is the time required to achieve a transport 
fidelity of at least 99% via resonant balancing. The black line 
(diamonds) is the time at which optimal fidelity is achieved 
for on-resonant transport without balancing. The fidelity is 
however quite low (~ 15%) in this case. The blue dashed line 
shows the number of PI centers in the network, N — 2. 



tween two separated NV centers via a bulk network con- 
sisting of randomly positioned nitrogen impurities (PI 
centers). Some experimental challenges remain: the lim- 
itation of pure state transfer requires that the system is 
at low temperature or that PI centers are indirectly po- 
larized by the NV centers; the requirement of using the 
isotropic XY Hamiltonian requires additional control for 
it to be generated from the natural dipolar Hamiltonian 
(for example via a combination of gradient and TOCSY 
mixing [21, 36, 39]). Still, considering that dephasing 
times in excess of 100/xs are routinely achievable [37, 38] , 
the balanced transport scheme may be experimentally vi- 
able for quantum communication in these architectures. 

Conclusion - In this paper, we showed that perfect 
quantum state transfer can be engineered even in the case 
of arbitrarily complicated network topologies, if the ends 
of the network are weakly coupled to the bulk. Transport 
speed can be improved by bringing the end spins on res- 
onance with a common mode of the bulk network. Alter- 
natively, it is possible to achieve unit transport fidelity, 
with lower control requirements, but on a longer time 
scale, by detuning the end-spins off-resonance. These 
transport strategies may allow the interlinking of quan- 
tum registers in a quantum information processor with 
very relaxed fabrication requirements. 
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Appendix A: Transport for end-spins on-resonance with a degenerate mode 



In the main text we described transport when the end spins are on resonance with an eigenmode of the bulk with 
eigenvalue EJ^ . Here we generalize to the situation where the eigenvalue is degenerate, with eigenvectors \am)- The 
projector onto the degenerate eigenspace is then, 



P = (|1) (1| + \N) (ATI) + J2 ("ml = (|1) (1| + \N) {N\) + Pm, 



(Al) 



meM 



and, to first order, transport is driven by the projection of the adjacency matrix into the this subspace, i.e., 
j^d _ pt^ep Now we consider that and A^ have the following forms, where "x" denotes a non-zero element: 



A'' = 



X X X X 

X X 

X X 

X X 

X X 

X X X X 



XXX 

X XX 
XX X 
XXX 



It follows that 



yielding 



= ; {N\A'\N) = ; {am'\A'\am) = 0, 



A' = ^(l|A<=|a„)(|l) {am\ + (1|) + ^(iV|A^|a„)(|iV) (a„| + |a„) (iV|), 



(A2) 



(A3) 



(A4) 



where we have used the fact that A^ is real, and hence (IjA'^lam) = (am|A^|l). Let us define the end connection 
vectors. 



Then, 



A^|l) = |ni) ; A'\N) = \nN). 
= + \am) (1|) + Y,{nN\am){\N) {ami + |am) {N\) 

m m 

= 5^^(ni|a„) ((a„|i) |1) 01 + {j\am) \j) (1|) + ^^(n;v|a„) ((a„|i) \N) {j\ + {j\am) |i) (A^|) 



(A5) 



3 m 



or simplifying. 



A'^ = J2iS,j\l){j\+SjM\N) {j\+h.c.), 



(A6) 



where (5(i_Ar)j = (?^(i.7v)l-PM|i) is the overlap of the end vector ni and the node j in the resonant subspace. 

To achieve balanced on-resonance transport we require that dij = Sjn for all j, which implies that both the 
end-vectors have equal projections in the resonant subspace. 



Pm \ni) = Pm Itin} ■ 



(A7) 



Appendix B: Transport Fidelity for A-networks 

Here we derive the maximum transport fidelity for a A-type network. A-type networks are interesting because the 
effective Hamiltonian of more complex networks reduces to A-network Hamiltonian in the weak-coupling regime (as 
shown in the prcvioiis section), but we will consider here the general case. 

For any network of adjacency matrix A, the fidelity function F(t) = \{N\ e~*"** has a simple expression in the 




FIG. 6. (Color online) A general detuned A-network with multiple A-paths between the end spins. Each leg of any of these 
A-paths may have an arbitrary detuning. For the representative case of the A-path going through node j (red), these detunings 
are Sij and SjN- 



eigenbasis of A, 



Fit) = ^(^1 A^) {N\k) {k\l) {l\i) cos [Ek - Ee) t. 

k,i 



This shows that the fidehty can be written as the sum 

Fit) ^ ^ Wn cos(/„t), with /„ = Ek- Ee, 



(Bl) 



(B2) 



that is, the frequencies are differences between eigenvalues of A, for which the corresponding eigenvectors \£) and \k) 
have non-zero overlap with |1) , |A^). 

We now consider a general A-network with multiple A-paths that connect the end-spins (see Fig. 6) and we will 
restrict the analysis to the adjacency matrix obtained in the on-resonance case only later. We write the adjacency 
matrix in terms of the coupling strength Sij and Sjn between the end spin and each j*'* spin in the bulk, which form 
the A path: 



A'^^Y. ' - (l^'Xll + + ^3Ni\j){N\ + \N){j\). 



(B3) 



Our strategy for finding the transport fidelity in A-networks is to first determine the eigenvalues of the adjacency 
matrix in Eq. (B3) and hence the possible frequencies at which information transport can occur. Then, we will use a 
series expansion to find an explicit expression for the fidelity. 



1. Frequency of Transport 



The eigenvalues of the adjacency matrix in Eq. (B3) are 

Ao = 0, {N — 4) degenerate 



where we defined, 



A3,4 - ±\/s^Ws^^, 

j 3 
j<k j<k 

6^ = YSij6jN 



(B4) 



(B5) 



While for a general network the frequencies of transport are differences between the eigenvalues of A"^, here there 



8 



are only four distinct frequencies because of the symmetries in the eigenvahies: 

/o = 0, /i=2Ai, /2 = 2A3, 



(B6) 

/3 = Ai - A3 = V2(^'-A2), /4 = Ai + A3 = V2(^' + A2). 



2. Series expansion 



With the frequencies found above, equation Eq. (B2) reduces to F{t) — X]i=o '^'^^ (/«^)- ^^'^ parameters 
Wi, we equate the Taylor expansion of Eq. (B2) and of the fidelity F{t) — \{N\ e~*"^* We only need the first five 

even power coefficients to fully determine {wi}, giving the series of equations, 

^«;, = |(7V|l)p = 

i=0 

4 

i=l 

i,E^^/'=4l(^l(^')'|l)l'=^4 (B7) 

1=0 

^ E^Jf = -^Re[(iV| {A-'Y \1){N\ {A'')' \l)]^C, 

i=Q 

^E"^^/' = ^!l(^l(^')'|l>l' + ^I^^KA^I 

■ i=0 




FIG. 7. (Color online) Positive frequencies of transport fidelity for a A-network (Fig. 6) consisting of two A-paths, 1 ^ j ^ N 
and 1 — > fc — > A'^. The frequencies are plotted as a function of the relative detuning between the paths, A^j, = {SijSkN — SjNSik)- 
The actual transport also contains symmetric negative frequencies and a DC (zero frequency) component. In general there are 
four frequencies of transport, derived in Eq. (B6). Note that when 5jM /5-ij is a constant for both paths, there are only two 
frequencies, S and 25, that carry the transport. 
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The expectation values can be evaluated exactly, yielding 



For general frequencies fi, the coefficients Wi are 

4 



2:4 



12 



-S^S 



(B8) 



720 



52 [95'^ - A2] 



Wo 



E 



Using the expressions for the frequencies in Eq. (B6), we find their explicit expressions in terms of 5, A and 8: 



Wo 



= W4 = -Wo 



(B9) 



(BIO) 



The fidelity is thus further simplified to 



m = gTZTAi h (V(5'+A2)/2) sin (^(^^ - A2)/2 



(Bll) 



3. Fidelity for random and degenerate networks 



Consider the case when the number of nodes is large, and the detunings 6ij and SjN are sampled from the same 
distribution, as it would be in a random network. Then, 



since the second moments of the random distribution should be equal. In this situation we have 

2 



- (A^ + 5^) = 



1 



E(^?.-'5f-) 







Then the condition 



5'' = S^- A^ 



is satisfied and the fidelity becomes 



(B12) 



(B13) 



(B14) 



F{t) = [sin + A2)/2j sin [ty/ {S^ - /2j\^ (B15) 

In the case of resonance to a non-degenerate mode, we have A = and the fidelity can reach its maximum F{t) — 
sin{St/V2Y = 1, for t = 'K/y/2S. 

For the case of interest in this work, a network where the end-spins are on resonance with a non-degenerate mode, 
the adjacency matrix of relevance in the weak regime is the reduced adjacency matrix, A'^. As shown above, in this 
case we have A = and the fidelity reduces to 



s 



St_ 

71 



n 4 



(B16) 



thus maximum fidelity can be reached only if the condition Eq. (B14) is satisfied. For example, the mirror-symmetric 
case 5ij = SjN, Vj yields the optimal fidelity F = I since then 5 = S. 
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Appendix C: Estimating matrix norms for different network topologies 
1. Different kinds of networks 



In this section, we consider different classes of networks, and estimate the norms of the corresponding adjacency 
matrices A. As described in the main text, the matrix norm of the bulk adjacency matrix is important in predicting 
the transport time. For example, the scaling of transport time is linear with 7 in the on-rcsonancc case; and the value 
of 7 implicitly depends on the norm of the bulk matrix. Hence a large bulk matrix causes an intrinsically high 7, 
and reduces the control requirements on the end-spins. All the networks considered are of N spins, and hence the 
adjacency matrices are N x N matrices. 

1. Random network: The matrix A has random entries in the range [0, 1] (with appropriate symmctrization). 
The random entries follow a uniform distribution, with no site-to-site correlation. Overall, this case represents 
a rather unphysical scenario, but will be useful in the computations that follow. 

2. Random network with scaling: A contains random entries from a uniform distribution scaled by 

l/{hd)^, where h is the Hamming distance between two nodes. It represents a network similar to a spin chain 
where all neighbor connectivities are allowed, and there is a possible spread in the position of the nodes from 
their lattice sites. 

3. Dipolar scaled regular (symmetric) network: We consider the network to be regular (symmetric) in two 

and three dimensions. With an appropriate choice of basis, this can be converted to a Bravais lattice. Special 
cases of interest arc the graphenc (honeycomb) lattice and the CNT (rolled honeycomb) lattices. 

4. Dipolar scaled regular network with vacancies: Here we consider the regular network above and introduce 
vacancies that are binomially distributed with parameter p. This approximately maps to the NV diamond 
system, where we consider transport through a PI lattice. 



2. Mathematical preliminaries 

1. The generalized adjacency matrix A of a network consists of positive numbers in the range [0, 1]. The matrix is 
symmetric and Hermitian. 

2. For the norm, we will use the Frobenius norm, which is the generalized Euclidean norm for matrices. 



N 
i 

where ai are the singular values of A. 



\ 



N 

El 



3. Random network 



Consider the right triangular form (R-form) of A, 



X X X X X 

X X X X 
X X X 

X X 
X 




(C2) 



Here x refers to random numbers uniformly distributed in the range [0, 1]. We have = Var[X] -|- [ElX])"^ = 

1/12 + 1/4 = 1/3. The total number of elements in the R matrix is. 



N 



N{N -1) 



(C3) 
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Hence, since the random numbers are assumed to be uncorrelated from site to site, we have, 

im - ^1^^ (C4) 

Fig. 8(a) shows the linear scaling of the norm in Eq. (C4), compared to the numerically obtained average of 100 
manifestations of random networks. The highest eigenvalue of A, E-^n^^ also scales linearly with N . 



4. Dipolar coupled random network 

Now we consider the case where there is 1/d^ scaling with the Hamming distance between two nodes. This represents 
a network similar to a spin chain where all neighbor connectivities are allowed, and there is a possible spread in the 
position of the nodes from their lattice sites. The adjacency matrix has the form, 





{2dp (3£i)3 (4d)3 (5d)3 

X X X X 

^ S3 (25p" (33p" (lap" 

" d3 {2d)3 (3d)3 

— ^ 

^ d3 (2d)3 

— 




(C5) 



As before, assuming that the sites are uncorrelated for the uniform distribution of random numbers x , we have, 



\Af = 



N-1 N-2 



1 



d6 (2d)6 (3d)6 

i=i ■' J=l 



Consider that. 



N-l 



TT 

945 ~ 

N-l 



1.01734 



1.036 



and the convergence is very rapid, i.e. it is true even for small N. Then, 



(C6) 
(C7) 



(C8) 



(C9) 



Fig. 8(b) shows that the scaling matches very well with the numerically obtained average norm of 100 mani- 
festations of random dipolar networks. The highest eigenvalue i?max approaches a constant 1.6/(P. 



5. Dipolar coupled regular network 

Here we consider a regular, symmetric network in two or three dimensions. To a good approximation, we can 
assume. 



\A\\ = n\\A\\cen 



(CIO) 
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(a) Random network 



(b) Dipolar network 




^0 20 40 60 80 100 "-"0 
Size of network N 




th 



E, 
E 



th 



50 100 150 200 
Size of network N 



FIG. 8. (Color online) Figure shows the scaling with network size A'^ of the matrix norms and largest eigenmodes of the 
adjacency matrices corresponding to (a) a random network and (b) a dipolar random network. The solid lines are average 
values obtained from 100 manifestations of the networks. The dashed lines are theoretical results. For the dipolar network, the 
largest eigenmode i?max approaches a constant LQ/cP (dashed magenta line). 

where ||y4||coU is the adjacency matrix of the unit cell of the underlying lattice, and n is the number of tihngs of this 
unit cell, 



N 



(Cll) 



where A^coii is the number of nodes per unit cell. 

Il^llccii depends on the choice of lattice in the particular network. Let us consider the case of a honeycomb lattice, 
where we assume only nearest neighbor interactions. This network is found naturally in graphene and CNTs. Then, 
1 1 A| I cell = 24/d'^. Hence, for graphene, 



1^11 = 



2^N 



(C12) 



6. Dipolar coupled regular network with vacancies 



Let the probability of a vacancy occurring be p. Once again we assume a binomial distribution. We also assume, 
that we can estimate the norm in this case by using tiling - i.e. we consider the vacancies only in the unit cells. 
Consider for simplicity the special case of graphene. For j vacancies, we have. 



The corresponding adjacency matrix. 



Hence the mean, 



\\A\\, 



4iVccii - j 



4jVcell(l~p) 



(C13) 



(CM) 



(C15) 



Hence, 



l^ll 



2^N(l-p) 

d3 



(C16) 



